Sticky grains do not change the universality class of isotropic sandpiles 
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We revisit the sandpile model with "sticky" grains introduced by Mohanty and Dhar [Phys. Rev. 
Lett. 89, 104303 (2002)] whose scahng properties were claimed to be in the universality class of 
directed percolation for both isotropic and directed models. Simulations in the so-called fixed- 
energy ensemble show that this conclusion is not valid for isotropic sandpiles and that this model 
shares the same critical properties of other stochastic sandpiles, such as the Manna model. These 
results are strengthened by the analysis of the Langevin equations proposed by the same authors to 
account for this problem which we show to converge, upon coarse-graining, to the well-established 
set of Langevin equations for the Manna class. Therefore, the presence of a conservation law keeps 
isotropic sandpiles, with or without stickiness, away from the directed percolation class. 
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I. INTRODUCTION 

Toy models of sandpiles are the archetypical examples 
of self-organized criticality emerging out of time-scale 
separation 1113,13. Sandpile models come in many differ- 
ent flavors (deterministic or stochastic rules Q, discrete 
or continuous variables 0, with or without height re- 
strictions 1^], etc.), but they usually consist in adding 
grains one by one until a local threshold (typically a con- 
dition on some slope or height) is reached, triggering a 
series of redistribution events, i.e. "avalanches" , which 
may lead to dissipation of sandgrains at the open bound- 
aries. Their numerical study is notoriously difficult and 
first led to a largely unsatisfactory situation in which 
"microscopic details" were believed to influence scaling 
Toperties, in contradiction with universality principles 
I . Major progress in favor of universality came when 
sandpile criticality was put into the broader context of 
standard non-equilibrium absorbing-state phase transi- 
tions |lliS|il|i2 [11 El El El Usui. 

Indeed, switching off both dissipation (open bound- 
aries) and driving (slow addition of grains), the total 
amount of sand or "energy" is conserved, and becomes a 
control parameter for these "fixed energy sandpiles" . For 
large amounts of sand, the system is in an active phase 
with never-ending relaxation events, while for small en- 
ergies it gets trapped with certainty into some absorbing 
state where all dynamics ceases (all sites being below 
threshold) . Separating these two regime s there is a criti- 
cal energy which was shown [la. ll4L Il5| to coincide with 
the stationary energy-density in the corresponding orig- 
inal sandpile (with slow-driving and boundary dissipa- 
tion). In this way, the exponents characterizing sand- 
piles can be related to standard critical exponents in 
an absorbing-state phase transition |l8l | (an alternative 
route not discussed here is to map sandpiles into stan- 
dard pinning-depinning interfacial phase transitions) . 

Using this approach, it was determined that stochas- 



tic sandpiles '2^1 generally do not belong to the di- 
rected percolation (DP) class, prominent among absorb- 
ing phase transitions, but to the so-called "conserved- 
DP" or Manna class (hereafter C-DP/Manna) charac- 
terized by the coupling of activity to a static conserved 
field directly representing the conservation of sandgrains 
[H [11 [H [H [M [ll- The field theory or mesoscopic 
Langevin equations describing this class reads: 



dtp = ap — bp^ 
dtE = DeS/^P 



LupE 



(1) 



where p(x, t) is the activity field (characterizing the den- 
sity of grains above threshold), i?(x, i) is the locally- 
conserved energy field, a,b,uj,D,a, and De are param- 
eters and ?7(x, t) is a Gaussian white noise. Equation 
represents a robust and well-established universal- 
ity class including not only stochastic sandpiles, but also 
some reaction-diffusion systems [l3, [2l|, [23 

In a recent Letter however, Mohanty and Dhar [2^ 
have argued that generic sandpile models with "sticky 
grains" or "inertia" , where some grains remain stable 
even if the local threshold is passed, should be in the 
DP class. The authors presented convincing numerical 
and analytical evidence that indeed this is the case for 
a directed two-dimensional system, which happens to be 
mappable into an effective one-dimensional directed site- 
percolation dynamics. Also, for isotropic (undirected) 
models with stickiness it was claimed that DP scaling 
holds. For this second case the authors presented a 
Monte Carlo simulation and, additionally justified their 
findings by arguing that the "right" set of Langevin equa- 
tions for inertial sandpiles should include a coupling of 
the form ujpQ{E — p — Ec) where 6 is the Heaviside step 
function and Ec is the instability threshold, substituting 
the bilinear coupling ujpE in Eq. 

The logic behind such a term is, in principle, reason- 
able [ill and it is argued in [2^ that considering this 
coupling, one should leave the Manna universality class 
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and return to DP, i.e. the conservation law should be ir- 
relevant in the presence of "stickiness" (see the schematic 
diagram in figure 4 of H^). 

In this paper, we argue that this claim in |23j | is un- 
founded: even in the presence of inertia/stickiness, the 
generic universality class of (undirected) stochastic sand- 
piles models remains the C-DP/Manna one. The paper 
is organized as follows: we first report on extensive sim- 
ulations of the model studied in [23 in the fixed-energy 
ensemble, from which we conclude that isotropic sticky 
sandpiles are in the C-DP/Manna class. In section III, we 
integrate the set of Langevin equations with a function 
coupling, which we find to be also in the C-DP/Manna 
class, and we perform a numerical renormalization treat- 
ment to show that these Langevin equations evolve upon 
coarse-graining towards Eqs.lQJ. 
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II. MICROSCOPIC "STICKY" SANDPILE 
MODELS 

The model proposed in |2^ is a variation of the Manna 
model: a discrete sandpile, defined on a one-dimensional 
lattice, with a height threshold slow sand addition 
and dissipation at the boundaries, but including a (stick- 
ing) probability 1 — p for grains to remain stable even if 
they are above threshold. Here we consider only the limit 
which possesses a critical point, i.e. the bulk-dissipation 
rate is set to zero. Following the strategy in |lH HI. . 
we analyze the "fixed energy version" of the model: we 
suppress grain addition and boundary dissipation, fix the 
total energy to E, and use p as a control parameter. Note 
that, owing to the existence of a non- vanishing sticking 
probability, arbitrarily large heights are allowed. Active 
sites (at which h > h^) are updated in parallel with the 
toppling occurring with probability p. 

We have implemented two different versions in which 
each toppling event redistributes two grains to the two 
nearest-neighbor sites, either randomly (stochastic redis- 
tribution rule) or regularly, with one grain onto each 
neighbor (deterministic rule). The methodology followed 
is standard for absorbing phase transitions: first, the crit- 
ical point is determined by studying the decay of activity 
from some initial condition varying p in a large system: 
for large p values, activity saturates (active phase), while 
for small p activity vanishes (absorbing phase). At the 
critical point, separating these two phases, Pc, activity 
decays asymptotically as a power-law with the critical 
exponent 9 = l3/v\\. For the stochastic and the deter- 
ministic rules, we find, respectively, pc = 0.84937(2) with 
9 = 0.120(8) and pc = 0.76750(3) with 9 = 0.115(8) 
(Fig. n^jb) . These estimates of 9 are in good agreement 
with the best evaluations for the Manna class in one di- 
mension, i.e. 9 = 0.125(2) |23], and clearly incompatible 
with the DP value 9 w 0.159. 

Next, using the critical value determined above, the 
variation of the stationary saturation value of activity 
at the critical point for smaller system-sizes is recorded. 



FIG. 1: (a-b) Log-log plot of the time-decay of the order pa- 
rameter (activity density) for different system sizes (from top 
to bottom: L = 64,128,256,512,1024,2048 and L = 2^**) 
and for the (a) stochastic rule and (b) the deterministic one, 
at their corresponding critical points pc = 0.84937(2) and 
Pc = 0.76750(3). From the slopes, we determine 9 = 0.120(8) 
and 6 — 0.115(8) respectively (DP slopes are plotted for com- 
parison). In (c-d) we plot the saturation values at the pre- 
viously determined critical points for the stochastic (c) and 
deterministic (d) rules respectively. From the scaling at the 
critical point we determine — 0.22(1) for both of them. 



From the expected scaling law Pst{p — Pc) ^ L , we 
determine — — 0.22(1) for both the stochastic and the 
deterministic rule as shown in Fig.[T|;,d. Again, this value 
is in good agreement with available estimates for the C- 
DP/Manna class ^ = 0.215(5) [Mllil, and incompati- 
ble with the DP value ^ w 0.252. 

We have also performed spreading experiments [Tolll^ 
(not shown) by following the standard procedure: we 
perturb a natural absorbing state (one generated by the 
system dynamics) to generate a small amount of localized 
activity and analyze how it spreads out at the previously 
determined critical point. We measured 77 = 0.39(3), 
S = 0.167(5), and z = 1.39(3) exponents for the number 
of active sites, surviving probability and average square- 
radius critical, respectively 0, El [13. These values 
are in good agreement with the best estimations for the 
Manna class and differ from their corresponding DP 
values ( 77 « 0.313, S w 0.159, and z « 1.258 |^). 



III. NUMERICAL STUDY OF LANGEVIN 
EQUATIONS 

After revisiting microscopic sandpile models with 
sticky grains, their critical behavior appears to be fully 
compatible with those of the C-DP/Manna class. We 
now turn to a study of the coupled Langevin equations 
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proposed by Mohanty and Dhar to describe their coarse- 
grained dynamics. The stochastic equations proposed in 
are Eqs. Q i.e. those of the C-DP/Manna class, 
except for the couphng term ujpE which is replaced by 
uipQ{E — p~ Ec) where Ec is the (microscopic) toppling 
threshold. The presence of this microscopic feature and 
of the step function O is surprising in so far as Langevin 
equations are usually understood as resulting from some 
coarse-graining of microscopic dynamics. In particular, 
the step function is unlikely to be a robust mesoscopic 
description, as it will be modified (probably transformed 
into a smoother function) upon coarse-graining. 

Discontinuous functions are notoriously difficult in the 
framework of renormalization group analysis. Moreover, 
even the "simple" equations resist standard pertur- 
bative renormalization attempts . The only available 
strategy then is direct numerical integration. The pres- 
ence of the (square-root) multiplicative noise term makes 
this a priori difhcult , but this technical difficulty was 
recently circumvented by the fast and quasi-exact sam- 
pling method introduced in [23|. 

We have used this scheme to numerically integrate the 
equations proposed by Mohanty and Dhar. These simu- 
lations yielded two sets of results: following the protocol 
recalled above, we studied the absorbing phase transition 
observed when varying the linear coefficient. We also in- 
troduced a local effective "mass" coefficient 

acff(x, t)=a + uj QiE{x, t) - p(x, t) - E,), (2) 

and studied its behavior upon coarse-graining in numer- 
ical simulations to analyze the relation between the two 
sets of Langevin equations. 



Phase transition 




b) 



FIG. 2; Direct numerical integration of the Langevin equa- 
tions proposed in |23| . (a) Decay experiments at the critical 
point a = 0.72308 (other parameter have been fixed as: sys- 
tem size L = 2^°, {E{x,t = 0)> = 0.5, 6 = 1, = 0.5, 
D = De = 0.25, Lo'^ = a'^ = 2, integrated with time-step 
dt — 0.25). (b) Finite size scaling at criticality (more details 
in the text). 
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Starting from an homogeneous, active, initial condi- 
tion, we studied the time-decay of spatially-averaged ac- 
tivity varying the control parameter a. As expected, alge- 
braic decay is found at the critical value separating expo- 
nential decay (absorbing phase) from saturation (active 
phase). The estimated decay exponent 6 = 0.130(5) is 
in perfect agreement with the C-DP/Manna class value 
(Fig. |2t). At the critical point, a — 0.72308, the scaling 
of the stationary activity for finite size systems yields the 
estimate (3/i'± — 0.22(1), (Fig. I^b) again in agreement 
with the C-DP/Manna value. These estimates are thus 
incompatible with the DP class values. Also, spreading 
experiments (not shown) fully confirm this result. 

B. Numerical coarse-graining 

The above results are easily understood when observ- 
ing the behavior of the Mohanty-Dhar Langevin equa- 
tions coarse-grained numerically. To do this, we build 
scatter plots of (ocfi) at vs the field difference {E — p)^, 



FIG. 3: Effective mass as defined by Eq. |2j as a function of 
the field difference averaged in Kadanoff blocks of size A'^, for 
a = 0.72313 in the active phase and lu = \/2. The vertical line 
corresponds to the threshold value he = 0.9. Observe that the 
larger the block size, the smoother the effective- mass dependence 
on the coarse-grained field difference. 



where the averages are taken on Kadanoff-blocks of 
length N. For = 1 (Fig.|21i), we obviously observe the 
0-function form. For A^ = 4 (Fig.lJb), the effective cou- 
pling can take intermediate discrete values between a and 
a+Lo, depending on the number of above-threshold micro- 
scopic sites in the block. When the block size is larger and 
larger (Fig. |3|:,d), a smooth function gradually appears. 
By retaining just the two leading terms in a Taylor ex- 
pansion of such an analytical function around the origin, 
we recover, at a coarse-grained level, the Langevin equa- 
tion for the Manna class Eqs.lQJ, i.e. a linear coupling 
term and a correction to the linear term, a in Eqs.(^. 
Higher order terms in the Taylor expansion can be ar- 
gued to be irrelevant from standard naive power count- 
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ing arguments. Therefore, it is not surprising that the 
set of Langevin equations including a Heaviside O func- 
tion should exhibit the same asymptotic behavior as the 
original C-DP/Manna class Langevin equations, Eqs.l^l. 

IV. CONCLUSION AND DISCUSSION 

To sum up, introducing "stickiness" in isotropic sand- 
pile models does not change their universality class, 
which remains that of the Manna model. We reached 
this conclusion via numerical simulations of microscopic 
models and Langevin equations proposed in [2^. We 
showed in addition that these Langevin equations "flow" 
towards those of the Manna class under some numerical 
coarse-graining procedure. 

One can wonder what is the reason why this conclu- 
sion does not hold for the directed (or anisotropic) sand- 
piles studied also by Mohanty and Dhar in ^] (see also 



[2^), which they proved to be in the DP class. These 
directed models, defined on a two-dimensional lattice in- 
clude an isotropic direction and a fully anisotropic one, in 
the sense that sand goes "downwards" in that direction 
but not "upwards" . This makes it possible to map the 
problem on DP in (1 + 1) dimensions, i.e. the anisotropic 
dimension can be taken as "time" . The local conserva- 
tion of energy is present also in these models, but "local" 
here means in "space-time" neighborhoods, while energy 
is not conserved in the isotropic spatial direction. Hence, 
Manna behavior does not appear, while DP scaling, the 
usual one in the absence of spatial energy conservation, 
is expected to emerge, as indeed was proved in 
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